Global maps of lunar neutron fluxes from the 
LEND instrument 


1. M. L. Litvak 1 ’*, 

2. I. G. Mitrofanov 1 , 

3. A. Sanin 1 , 

4. A. Malakhov 1 , 

5. W. V. Boynton 2 , 

6. G. Chin 3 , 

7. G. Droege 2 , 

8. L. G. Evans 4 , 

9. J. Garvin 3 , 

10. D. V. Golovin 1 , 

11. K. Harshman", 

12. T. P. McClanahan 3 , 

13. M. I. Mokrousov 1 , 

14. E. Mazarico 3 , 

15. G. Milikh 5 , 

16. G. Neumann , 

17. R. Sagdeev 5 , 

18. D. E. Smith 3 , 

19. R. Starr 6 , 

20. M. T. Zuber 7 

Article first published online: 5 JUN 2012 

DOI: 10. 1029/201 1JE003949 

Copyright 2012 by the American Geophysical Union 

Issue 



Journal of Geophysical Research: Planets (1991-2012) 


Volume 1 17. Issue El 2. December 2012 


Additional Information fShow All) 


How to CiteAuthor InformationPublication History 

SEARCH 

Search Scope 


Search String 

. AGU EASI Search > 

• Browse by Index Term > 

ARTICLE TOOLS 

. Get PDF (748 IK) 

• Save to My Profile 

• E-mail Link to this Article 

• Export Citation for this Article 

• Get Citation Alerts 

• Request Permissions 

More Sharing ServicesShare | Share on citeulikeShare on connoteaShare on deliciousShare on 
www.mendeley.comShare on twitter 

• Abstract 

• Article 

• References 

• Supporting Information 

• Cited By 

Get PDF (748 IK) 

• Go here for SFX 


Keywords: 

. LEND; 

• Moon; 

• epithennal neutrons; 

• hydrogen; 

• water 


Abstract 


[1] The latest neutron spectrometer measurements with the Lunar Exploration Neutron Detector 
(LEND) onboard the Lunar Reconnaissance Orbiter (LRO) are presented. It covers more than 

1 year of mapping phase starting on 15 September 2009. In our analyses we have created global 
maps showing regional variations in the flux of thermal (energy range < 0.015 eV) and fast 
neutrons (>0.5 MeV), and compared these fluxes to variances in soil elemental composition, and 
with previous results obtained by the Lunar Prospector Neutron Spectrometer (LPNS). We also 
processed data from LEND collimated detectors and derived a value for the collimated signal of 
epithermal neutrons based on the comparative analysis with the LEND omnidirectional detectors. 
Finally, we have compared our final (after the data reduction) global epithermal neutron map 
with LPNS data. 

1. Introduction 

[2] The Moon is continuously bombarded by high energy charged particles (primarily protons 
and alpha particles) from Galactic Cosmic Rays (GCR). The result of such interactions is 
production of neutrons over a wide energy range within the upper (~ 1 m) subsurface. The 
observation of the neutron leakage spectrum may provide information about the hydrogen 
content of the Moon itself, as well as about the space environment (for example variations of 
GCRs, see Mitrofanov et al. [20091) . 

[3] Regional variations of neutron flux in the thennal energy range (neutrons with energies 
below 0.4 eV) correlate with concentrations of major and minor lunar soil forming elements (like 
Fe, Ti, Gd, Sm) having large macroscopic absorption cross sections. For example, the 
distribution of Fe and Ti is highly “non-homogeneous” between the nearside mare basins and 
farside highlands and as a result produces the significant regional variations of thermal neutrons 
across the lunar surface [see, e.g., Lawrence et a /., 2002 ; Elphic et al., 1998, 2000, 20021 . The 
epithennal energy range (neutrons with energies from 0.4 eV up to 100 keV) is another key 
subsurface composition marker because the intensity of neutron flux in this energy range is 
sensitive to the abundance of hydrogen. Even a small amount ( — 1 00 ppm) of H causes 
significant depression of epithermal neutron flux. Finally, fast neutrons (neutrons with energies 
from 100 keV up to 15 MeV) may be used to analyze the average composition of lunar soil via 
correlation of intensity of fast neutron flux and average atomic mass of the lunar soil. 

[4] More than 10 years ago lunar neutron leakage spectra were acquired by the Neutron 
Spectrometer instrument onboard the Lunar Prospector mission (LPNS; for instrument details, 
see Feldman et al. [20041) . LPNS was able to create global maps of lunar neutron albedo in 
different energy ranges including observations of thermal, epithennal and fast neutrons [ Feldman 
eta/., 1998a , 1998b ; Maurice et al., 20041 . The most intriguing result was the observation of 
regional variations of epithennal neutron flux across the lunar surface. LPNS revealed vast polar 
regions of epithermal neutron suppression (4-5% lower in comparison with surrounding low 
latitude areas with similar composition except H) at both the north and south poles, interpreted as 
local water ice rich areas (presumably within permanently shadowed regions) or as a reservoir of 
volatiles created by solar wind migrate onto cold poles [ Feldman et al., 1998b ; Crider and 
Vondrak, 2000 ; Feldman et al., 2000 , 2001 ; Lawrence et al, 20061 . The best spatial resolution of 


LPNS (full width at half maximum resolution —46 km at low altitude about 30 km [see Maurice 
et al., 20041 ) was too poor to resolve local areas with highest H abundance leaving open the 
question of whether the hydrogen was uniformly distributed within the extended suppression 
regions or localized in permanently shadowed regions (possibly from comet water ice 
deposition). 

[5] Global maps of thermal and fast neutrons also have shown significant regional variations of 
neutron flux across the surface closely related with non-homogeneity of soil composition. Thus, 
nearside mare basins and South Pole Aitken (SPA) basin are visible in global maps as strong 
depressions of thermal neutron flux [see, e.g., Feldman et al., 1998a ; Lawrence et al ., 20021 . The 
highlands on the farside of the Moon, on the contrary, are seen on a global map of thermal 
neutron flux as high intensity areas, poor in major absorbing elements like Fe and Ti \ Feldman et 
al. , 1998a1 . The significant difference in soil composition of nearside mare basins and farside 
highlands is also clearly visible in the fast neutron global map. It appears as a higher intensity of 
fast neutron flux in mare basins in comparison with lower intensity in highlands. The major 
cause of this difference in this case is a variation of average atomic mass of the lunar soil with 
larger values (more weight fraction of Fe) in mare basins and lower values (more weight fraction 
of Al) in highlands [see Gasnault et al 20011 . 

[6] The next mission to implement nuclear methods for studying lunar soil composition appeared 
in 2009 when another neutron spectrometer, LEND (Lunar Exploration Neutron Detector), was 
launched into lunar orbit onboard the Lunar Reconnaissance Orbiter (LRO) mission \Chin et al ., 
20071 . LEND consists of several neutron sensors measuring lunar neutron flux in different 
energy ranges with different spatial resolutions. It is similar to LPNS in the global monitoring of 
thermal, epithermal and fast neutrons by omnidirectional neutron detectors, which have spatial 
resolution strongly limited by the altitude, but LEND has the significant advantage of measuring 
collimated epithermal neutron fluxes with much better spatial resolution (see Table 1 and Figure 
1). LEND physical/geometrical resolution is significantly better than previous measurements 
made by LPNS. At the low altitude of —30 km LPNS claimed resolution of 46 km measured as 
full width at half maximum (FWFIM). During mapping phase LRO was flying at the altitude 

of —50 km. At this altitude LEND can provide FWFIM = 10 km. In December 2012 LRO was 
moved into the frozen elliptical orbit with periapsis about 30 km above the south pole. It gives 
possibility to map neutron flux by LEND collimated detectors with spatial resolution 
FWHM = 6 km. These estimations were verified through the numerical modeling and 
calibrations [see, e.g., Mitrofanov et al. , 2008 , 2010a , 2010b , 20111 . 

Table 1. The Description of LEND Neutron Detectors 
Detector Description 

Two ’He proportional counters which are attached to the opposite sides of bottom 
part of the LEND collimator body. They are oriented along the velocity vector of a 
spacecraft forming a Doppler pair of detectors. Each detector has efficiency in wide 
STN 1 and energy range equal approximately to 1 for thermal energies and decaying with energy 
STN2 for the epithermal neutrons. The neutron sensor facing forward the spacecraft 

velocity vector detects significantly more neutrons than detector facing in backward 
direction. This difference in counting rates strongly depends on neutron's energy 
having maximum at low thermal energies (<0.015 eV - equivalent energy for the 


Table 1. The Description of LEND Neutron Detectors 


Detector 


Description 


STN3 


SETN 


CSETNs 


SHEN-N 


neutrons with spacecraft velocity) and fast dropping down to zero for the epithermal 
energies (>0.4 eV). As result the difference in counting rates is given by thermal 
neutrons hitting the front side of spacecraft and free of any epithermal neutrons and 
thermal neutrons background produced in the spacecraft body. 

The omnidirectional 3 He proportional counter, which is attached to the upper part of 
collimator body. This sensor is identical to the Doppler proportional counters having 
the same sensitivity to the thermal and epithermal neutrons. This detector is used to 
support the estimation of spacecraft background produced by thermal neutrons. 

The omnidirectional 3 He proportional counter, which is attached to the upper part of 
collimator body (opposite side to the STN3). This sensor is covered with 1 mm sheet 
of Cd to stop registration of thermal neutrons with energies below 0.4 eV. The data 
from this detector are used to provide global mapping of epithermal neutrons with 
low spatial resolution (full width at half maximum is more than altitude of 
spacecraft). The difference of counting rates from STN3 and SETN characterizes the 
thermal component of lunar emission with the energy of neutrons below the Cd 
threshold at 0.4 eV. 

Four 3 He proportional counters, which are identical to STN1, STN2, STN3 and 
SETN. Tops of these detectors are open and directed to the nadir and sides are 
surrounded with polyethelene - 10 B collimator walls, bottom part is covered with 
cadmium and 10 B. It provides the angular efficiency function with half width at half 
maximum equal to 5 km at nominal 50 km altitude. These detectors are primarily 
used to create polar maps of the Moon with enhanced spatial resolution. 

This sensor is a stilben crystal scintillation detector with time separation circuit to 
split resulted detection signal into gammas and neutrons (registered by recoil protons) 
channels. This detector is integrated in the center of collimator and enclosed with 
anticoincidence plastic scintillator to remove «pollution» counting rate produced by 
space protons. It detects neutrons with energies from 0.5 MeV up to 15 MeV. The 
collimator walls are not efficient enough to stop fast neutrons as epithermal neutrons, 
but nevertheless it significantly depresses the flux of high energy neutrons leading to 
the narrowing field of view of SFIEN-N in comparison with omnidirectional 
detectors. The data from this sensor are used for global mapping of fast neutrons. 


Figure 1 . Instrument LEND has the set of nine sensors of thermal, epithermal and fast neutrons 
and collimation module. 



[7] The previous statement should be necessarily combined with the counting statistic measured 
by LEND. In other words, which portion of the measured counting rate came from the detection 


of epithermal neutrons in the open field of view (FOV) and which is induced by neutrons 
penetrating through the collimator walls or backscattered in the spacecraft body. It significantly 
influences the observational (exposure) time needed to distinguish small areas on the Moon 
surface at the statistically significant level of more than 3 er. Our analysis (see section 5) has 
shown that the total counting rate in the collimated detectors is about 5.1 counts per second. It 
consists of ~ 1.7 counts per second measured as a collimated signal and ~3.4 counts per second 
resulted from the spacecraft background and penetration through the collimator walls. In our 
paper, we only focus on the global neutron maps derived from the omnidirectional detectors and 
do not develop applications of this result to the study of local polar areas. But other publications 
Mitrofanov et al., 2010b , 20111 concluded that such proportion in the counting rate is sufficient 
to distinguish inhomogeneity in the distribution of neutron fluxes at small scale (large polar 
shadowed areas with 100-1000 km 2 ) at least during first year of mapping. 

[8] LEND detects neutrons starting from thennal energies (less than 0.015 eV) up to high energy 
(~15 MeV) neutrons. The main goals of this paper are to present the first global LEND maps of 
thermal, epithermal and fast neutrons corrected for various systematic effects; to study regional 
variations of neutrons in terms of soil elemental composition and to present results of a 
preliminary comparison with LPNS data. 

2. LEND Instrumentation, Data and Methods 

[9] LEND is a multidetector system consisting of 9 different neutron detectors responsible for 
registration of neutrons with different energies (for details, see also Mitrofanov et al. f2008, 
2010a]). All detectors can be divided into two groups. One is a group of detectors installed inside 
a heavy polyethylene - 10 B collimator. These detectors have a narrow field of view, which 
provides the ability to monitor lunar neutron fluxes of epithermal and fast neutrons with 
improved spatial resolution. Another group of detectors are omnidirectional sensors, which are 
located outside the collimator body and detect thermal and epithermal neutrons with a spatial 
resolution defined only by spacecraft altitude (>75 km FWFIM at the 50 km circular orbit during 
mapping phase). The short description of each LEND detector is presented in the Table 1 and 
Figure 1 . 

[10] In section 3 we present data derived from the two omnidirectional proportional counters 
(STN 1 and STN2) mounted on opposite sides of the collimator and aligned along the spacecraft 
velocity vector \ Mitrofanov et al ., 201 Oak Combined they act as Doppler pair detectors, sensitive 
to the registration of thermal neutrons in the energy range below 0.015 eV (see Feldman and 
Drake [19861 for detailed description of Doppler filter technique proposed for registration of 
thermal neutrons). Based on these data, we can create a global map of the thermal neutron flux to 
compare with the LPNS results. 

[11] In section 4 we present data derived from the fast neutron detector, which is an organic 
(stilben crystal) scintillator (SFIEN-N detector) with risetime discrimination that splits the signal 
into gamma and neutron (registered by recoil protons) channels [ Mitrofanov et al ., 2008, 2010a], 
This detector is located in the center of the collimator and surrounded with an anticoincidence 
plastic scintillator to remove the background counting rate produced by space protons. The 
whole detector assembly is inherited from the fast neutron detection technique used in the High 


Energy Neutron Detector (HEND), part of the Mars Odyssey instrument suite \ Boynton et al., 
20041 . Using this data set we have created a global map of fast neutron flux for comparison with 
the map provided by LPNS. 

[12] In section 5 we concentrate on the comparison of different data sets from omnidirectional 
SETN detector and four collimated CSETN1-4 detectors with the narrow field of view 

[ Mitrofanov et al . , 201 0a] . Our goal here is to derive and validate the epithermal component of 
neutron flux, create global maps and compare with the LPNS map. 

[13] All data discussed in section 3—5 are taken from the LEND official Planetary Data System 
(PDS) data for a time period from 15 September 2009 to 15 December 2010, a little longer than 
the official Exploration Systems Mission Directorate primary mission. 

[14] All data sets are passed through the LEND standard data reduction procedures that take into 
account various systematic effects including correction for efficiency changes, trend of GCRs, 
exclusion of strong Solar Particle Events, separation between charged particles and neutrons, 
altitude and temperature correction, and analysis of cruise data to identify and eliminate 
spacecraft background (see also M. L. Litvak et al., LEND neutron data processing for mapping 
of the Moon, submitted to Journal of Geophysical Research , 2012). 

3. Distribution of Thermal Neutron Fluxes 

[15] In this section we discuss the main results of global mapping of thermal neutron fluxes 
measured by LEND Doppler detectors. As it was already mentioned in section 2 , for this analysis 
we have used data from two omnidirectional proportional counters aligned along the spacecraft 
velocity vector and separated by the LEND collimator. The difference in the counting rates 
between these two detectors is proportional to the thermal neutron flux at the lowest energies 
(primarily the neutrons with energies below 0.015 eV). More than one year of primary mapping 
data (circular orbit with an average altitude of 50 km) was used as input for mapping the thermal 
neutron flux data. This time period was divided into four shorter time intervals delineated by 
changing the spacecraft orientation to the velocity vector. The mapping phase began with 
velocity vector oriented along the spacecraft's —x axis (The x axis can be described as a vector 
connecting STN2 and STN1 detectors facing in the direction of STN1, see Figure 1) . After just 
two weeks the spacecraft orientation was changed to the opposite one, with velocity vector along 
the +x axis, where it remained until the end of March 2010, then flipped back to the —x axis. The 
fourth and final period, with the velocity vector aligned with spacecraft's +xaxis, was started at 
the end of September and lasted up to the end of analyzed period of the mapping. The counting 
rates in two LEND Doppler detectors, STN1 and STN2, are “velocity direction” dependent with 
the lead detector measuring thermal and epithermal neutrons and the trailing detector measuring 
mainly epithermals. The thermal flux for each of the four time periods is the difference of the 
detectors, STN1-STN2 for the +xorientation and STN2-STN1 for the -xorientation. The data 
reduction for each Doppler detector included corrections for efficiency and variations in GCR 
flux. The time profiles of counting rates measured by both omnidirectional detectors SETN and 
STN3 and all collimated detectors show very similar long-term variations (a slow reduction of 
counting rates by ~15%) during the first year of mapping. It correlates well with the changes of 
GCR flux observed by GOES and ACE missions. STN1 and STN2 detectors also show the same 


trend of counting rates but are significantly modified by the numerous changes in spacecraft 
orientation. For detectors STN1, SETN, CSETNs and SHEN-N we have used the time profiles 
averaged for a long period of time (two weeks) to get corrections factors for the GCR variations 
(see Litvak et ah, submitted manuscript, 2012). For STN1 and STN2 we have used relative 
trends derived from the averaged STN3 time profile to make correction for the GCR long-term 
behavior. 

[16] Another correction was applied to account for the non-homogeneous distribution of 
spacecraft mass around the Doppler detectors. The field of view of STN 1 is partially obscured by 
the spacecraft body while the field of view of STN2 is clear of it. This means that some 
additional background Cb g d/smi (thermal neutrons from spacecraft) is added to the counting rate of 
STN 1 in comparison with the counting rate of STN2. When the velocity vector is oriented 
along +xaxis, STN 1 is the lead detector and STN2 trails. So the difference in counting rates, 
STN1-STN2, for this orientation is larger than difference in STN2-STN1, as determined by the 
spacecraft orientation in the —x axis direction. This is because in the first case, the spacecraft 
background is added to the difference of counting rates between STN 1 and STN2, but in second 
case it is subtracted making it smaller. Roughly, the background counting rate may be estimated 
from the comparison of two spacecraft orientations: 2 x Cb g d/stni = [(C st ni-C stn 2) - (C stn 2-C s tni)], 
where C s tni and C stn 2 are measured counting rates in the detectors STN1 and STN2. Taking this 
equation into account one may adjust all spacecraft orientation time periods to the standard 
counting rate in Doppler detectors subtracting Cb g d/ S tni from C stn i. Following this procedure we 
estimated counting rate in Doppler pairs for each time interval, and mapped it with resolution 1° 
x 1°. The resulting statistical uncertainty at each 1° x 1° pixel depends on the LRO exposure 
time, which is a function of spacecraft orbiting and instrument on/off periods. The average value 
of the relative error in each pixel is about 7%. The final map of counting rate is smoothed with a 
Gaussian filter of FWFIM = 60 km, and presented in Figure 2 . It is comparable with the spatial 
resolution of an omnidirectional detector at an altitude of 50 km (FWFIM ~ 75 km). 

Figure 2. Smoothed (Gaussian filter with FWFIM = 60 km) map (in Robinson projection) 1° x 1° 
of thennal neutron flux from lunar surface measured as a difference in counting rate (counts per 
second) between pair of LEND Doppler detectors (STN1 and STN2). 



[17] The thermal neutron flux intensity shown on the LEND map has significant dynamic range 
with a factor 3.5-4 and is reasonably well correlated with the main patterns of lunar soil 
composition. The major effect is seen on the nearside of the Moon in mare basaltic terrains 
where there are substantial depressions of thermal neutron flux. In contrast, the highlands on the 
farside of the Moon show significantly higher intensities of thermal neutron flux. These regional 
variations have been observed in previous experiments and were analyzed in details using data 
from LPNS and Lunar Prospector Gamma-Ray Spectrometer (LPGRS) instruments. The 
correlation of lunar maria with considerable reductions in thermal neutron flux intensity has been 
interpreted as an enhanced abundance of such primary soil forming elements as Fe and Ti and 
minor and trace elements like Gd, and Sm that have large neutron absorption cross sections. 
Comparison of the distribution of thermal neutron flux detected by LPNS and 7.6 MeV gamma 


ray doublet (Fe neutron capture energy band) counting rates measured by LPGRS has shown that 
a key factor defining the global distribution of thennal neutrons is iron (FeO) abundance. FeO 
abundance has been found to vary from ~5% by weight in the highlands up to 22-23% weight in 
western mare regions with intennediate values 8-10% at South Pole Aitken basin [see, 
e.g ., Lawrence et al., 20021 . This pattern correlates very well with the distribution of thermal 
neutron fluxes measured by the LEND instrument. So, the direct correlation (pixel-to-pixel) 
between LPNS (taken as special PDS product prepared by LPNS team according to data 
reduction procedures described in Maurice et al. [20041) and LEND maps of thermal neutron flux 
shows very high linear correlation coefficient 0.98 [see, e.g., Press et al ., 19921 . It is illustrated 
in Figure 3 . 

Figure 3. The correlation between LEND thermal neutron Doppler detectors counting rate (y 
axis, counts per second) and LPNS counting rate from thermal neutron detector (x axis, counts 
per second) using a 2° x 2° regular mapping grid. The linear best fit is shown by dashed line. 



[18] Figure 4 summarizes the results of a comparative analysis of regional variations of thermal 
neutrons measured on global scale by the LPNS and LEND instruments. 

Ligure 4. Orbital phase curves (see section 3) showing regional latitude variations of thermal 
neutron flux across the lunar surface. The x axis is the phase angle measured in degrees of arc 
from 0° up 360°, where 90° corresponds to south pole and 270° to the north pole. They axis is 
the normalized counting rate of fast neutrons measured by different instruments. LEND data are 
shown by solid curve and LPNS data by dashed curve. 


i 
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[19] In addition to the global map presented in Figure 2 we can characterize regional variations 
of neutron flux as a function of latitude. We have combined all orbital measurements into one 
average orbit, keeping latitude dependence and averaging different longitudes. Such orbital 
phase profile (OPP) presents a function of phase (latitude) angle (p measured in degree of arc 
from 0° to 360° (binning of full orbit around Moon). This is similar to averaging an ensemble of 
counting rates within wide polar bands, as suggested by Feldman et al. [200 1] . In this approach 
we derive one-dimensional view of the pattern of regional variations presented on the global map 
(Figure 2) by averaging counting rate across longitudes but leaving its dependence from the 
latitude. It partially suppresses (large areas like mare basalt terrains still will be quite visible) 
longitude dependence but significantly increase statistic of counts. Especially it is helpful to 
search for possible significant variations of neutron flux inside polar regions in comparison with 
low latitude areas. Our definition of average orbit can be described as follows: The range of 



phase angles (p from 0° to 90° corresponds to the averaging parts of all available orbits when the 
spacecraft was moving from equator to south pole within [-180°, 0°] longitudes. The range from 
90° to 270° correspond to averaging all orbits when spacecraft was moving from south pole to 
north pole within longitudes [0°, +180°] and finally, <p ranges from 270° to 360° correspond to 
the averaging parts of orbits when spacecraft was moving from north pole to equator within 
longitudes [-180°, 0°]. 

[20] To compare regional variations of thermal neutron fluxes observed by the LPNS and LEND 
instruments, we have created OPPs and nonnalized each curve by dividing by the average 
counting rate. It was done to exclude differences in the absolute values of counting rates 
measured by different instruments and leave only relative variations in OPP's shape from one 
lunar region to another. The nonnalization was performed in accordance with expression 

Cnnrm = C tca i/< C rca i> , where C r eai is the absolute value of counting rate measured by one of the 
instruments and <C rea i> is its mean value calculated for the whole lunar surface. From Figure 4 it 
is seen that the normalized LEND and LPNS OPPs have similar variations. Discrepancies are 
most visible in the polar and southern intennediate latitude regions. Poleward of 70° both in the 
south and north, the LEND profile shows higher values (in comparison with values observed at 
equatorial regions) of thennal neutron flux than in the LPNS data. At the intermediate southern 
latitude belt (30S-75S) the thermal flux measured by LEND is lower than that from LPNS . 

[21] We suggest that a possible explanation of the polar differences lies in the latitude variation 
of the lunar subsurface temperature. LPNS measured thermal neutron fluxes in the energy range 
below the cadmium threshold (<0.4 eV [see Feldman et al ., 2004]) derived from the difference 
between counting rates from a bare He proportional counter and one covered with 0.63 mm of 
Cd. The primary method for detection of thennal neutron fluxes by LEND is based on the 
Doppler technique, measuring thermal neutrons with considerably less energies (mostly below 
0.015 eV) than those measured by LPNS. Little et al. [2003] and Lawrence et al. [20061 have 
concluded that thennal neutron fluxes measured by LPNS are temperature dependent and vary 
with soil type. It was verified with numerical simulations that thermal neutron fluxes in this 
energy range increase with increasing temperature of the soil, leading to the maximum variations 
between equator (-250 K) and polar latitudes (—100 K) of about 5%. Thermal neutrons detected 
by LEND belong to the low energy tail of the Maxwellian distribution and follow an opposite 
temperature trend. With increasing temperature the peak of the Maxwellian distribution is shifted 
to higher energies. It decreases the number of neutrons in the energy range below 0.015 eV 
(where Doppler detector efficiency is maximum). This leads to decreasing counting rates with 
increasing soil temperature. As a result, if the normalized LPNS and LEND counting rates fit 
each other at the equatorial regions (with high surface temperature) then at polar latitudes (with 
lower temperatures) we should expect discrepancies in their regional variations: larger amplitude 
of variations for LEND (increased by temperature effect) and smaller one for LPNS (decreased 
by temperature effect). Preliminary simulations of LEND data show that temperature variations 
may be as large as 5-8%. In this paper we have not included a quantitative correction of the 
thermal neutron counting rates (measured by Doppler detectors) for soil temperature. It is a 
second order effect, much smaller than global regional variations of the thermal neutron flux. 

One can compare possible temperature induced variations <10% with full scale of regional 
variations of 300-400% presented on the map in Figure 2 . As a result this effect does not 


significantly alter the main patterns of the global map of lunar thermal neutron flux (high pixel- 
to-pixel linear correlation coefficient between LPNS and LEND maps shown in Figure 3) . 

[22] For southern moderate latitudes (30S-75S) we hypothesize that differences between LEND 
and LPNS OPPs are due to variations of soil composition and different sensitivity of neutrons 
with energies <0.015 eV (LEND Doppler filter) and neutrons with energies <0.4 eV (LPNS 
omnidirectional detectors) to detect such soil variations. 

[23] To investigate differences between LEND and LPNS OPPs at equatorial and intermediate 
latitudes in more detail we have perfonned pixel-to-pixel comparison of normalized (to adjust 
different instrument sensitivities) LEND and LPNS counting rates. The result of the comparison 
is presented as a global map of a value D of pixel-to-pixel difference. The value Dy at each pixel 
{i,j } is defined as: 


I LENDu LPNS „ 

n _ | vM>-Tnm 

U v- s ' J 

where LENDy is the LEND counting rate at the pixel {i,j},<LEND> is the LEND average 
counting rate (calculated for the whole lunar surface), LPNSy is the LPNS counting rate at the 
pixel {i,j},<LPNS> is the LPNS average counting rate (calculated for the whole lunar surface) 
and finally e)y is the total statistical error of the numerator in equation (1) deduced from the 
statistical uncertainties of LEND and LPNS counting rates at the pixel {i,j}. 

[24] Conceptually the definition of Dyis close to the definition of signal-to-noise ratio (SNR) 
showing how the value of signal is statistically different from the background. In our case the 
difference between normalized LEND and LPNS counting rates (nominator of (1)) is used as a 
“signal” and the statistical error of this difference (denominator of (1)) is a measure of “signal” 
significance. The value of Dy presents a number of statistical deviations (sigmas) between LEND 
and LPNS data. If it exceeds 3 sigma level the difference between LEND and LPNS in the given 
pixel can be considered as significant. The resulted map of Dy is shown in Figure 5 . The yellow 
and red colors on this map correspond to the most significant differences between LEND and 
LPNS regional variations. It is seen that such “discrepancy” areas are distributed at the poles 
(already discussed above) and at the equatorial and moderate latitudes with enhanced content of 
Fe (mostly at South Pole Aitken basin). 

Figure 5. The map of comparison (number of statistical deviations) between LEND thermal 
neutron counting rate (pair of Doppler detector) and thermal neutron counts from LPNS detector. 
For details, see section 3, equation (1) . 



[25] We may conclude that the difference between LPNS and LEND regional variations of 
thermal neutron fluxes at mid-southern latitudes is produced by a different sensitivity of thermal 
neutrons with different energies due to the variations of soil composition. 


4. Distribution of Fast Neutron Fluxes 


[26] In this section we present results of global mapping of the fast neutron fluxes measured by 
LEND SHEN-N detector. Here we have used data from the organic scintillator stilben 
(C6H5CH = CHC6H5) located inside the collimator (for details, see Mitrofanov et al. [2008 , 
2010a]). This detector is used for detection of high energy neutrons by recoil protons in the 
stilben crystal lattice. The incident neutron loses its energy, E, through the collision of hydrogen 
nuclei leading to a nearly random and unifonn spectrum of protons in the energy range from 0 to 
E. The LEND collimator is not as effective in stopping fast neutrons as it is for the epithermal 
neutrons, but it nevertheless narrows the SHEN-N field of view. The real spatial resolution of the 
detector measured as full width at half maximum is about 40 km. SHEN-N data are measured as 
a 16 channel spectrum of counting rate induced by fast neutrons in the energy range from 

~500 keV up to 15 MeV. To improve statistics we combine all higher channels (#10-15) 
excluding low energy neutrons near the detector threshold (to exclude thresholds systematic 
effects and low statistic channels) and the highest energy integration channel (#16), which sums 
all counts from recoil protons with energies above 10 MeV (to exclude very energetic neutrons 
which are insensitive to the soil composition). 

[27] It is known [e.g., Feldman et al., 1998a] that fast neutron flux is an important component of 
the lunar neutron environment, is sensitive to the soil composition, and may be used to infer 
average elemental composition. Multiple publications have confirmed the linear relation between 
fast neutron flux and average atomic mass of soil composition, using both experimental data and 
numerical simulations [ Maurice et al 2000; Gasnault et al 200 1] , 


[28] We have produced the LEND global map of distribution of counting rate from the fast 
neutron detector with a resolution of 1° x 1° degree. Pixels on the map have varied statistical 
uncertainty depending on the LRO exposure time. The mean relative error is about 1 1%. As in 
case of thennal neutrons map we have smoothed fast neutrons map by Gaussian filter with 
FWHM equal to 60 km (SHEN-N has better spatial resolution but significantly less counting 
statistics). The result is presented in Figure 6 . It clearly shows a non-homogeneous distribution 
of fast neutron flux across the lunar surface. The basaltic mare terrains at the low latitudes are 
seen as bright areas with substantially higher neutron flux than those observed in the highlands. 
The same pattern of fast neutrons is also seen in the LPNS data [ Feldman eta/., 1998a] , The 
observed effect is explained as a correlation of fast neutron flux with average atomic mass: in 
mare regions enriched with Fe and Ti the average atomic mass of the soils is larger than in the 
highlands (Al and Ca - dominating elements). To confirm this conclusion we have compared the 
average atomic mass calculated as equation (2) for kn own elemental compositions of lunar 
landing sites for Apollo/Luna missions and LEND counting rate of fast neutron detector. 





( 2 ) 


where A n its atomic mass of element n (in atomic mass units, amu), and w n its weight fraction. In 
Figure 7 one can see the results of such analysis. The x axis on this figure corresponds to the 
LEND counting rate derived for each landing site (smoothed counting rate in pixel 1° x 1° 
centered on a landing site). They axis corresponds to the average atomic mass calculated by 


equation (2) for each landing site. A combination of Apollo and Luna mission landing sites 
compositions was used [see, e.g., Vinogradov et al., 1966; Metzger et al . , 1973 ; Lawrence et al.. 
20061 . All eight data points (Apollo 11, Apollo 12, Apollo 15, Apollo 16, Apollo 17, Luna 16, 
Luna 20, Luna 24; the same set of data as analyzed in Gasnault et al. [20011) are well described 
by a simple linear fit (see the solid curve in Figure 7) . The linear correlation coefficient \ Press et 
al. , 19921 is considerably high, estimated to be ~0.9. The possible deviations from the linear law 
are explained by large statistical errors of the counting rate, poor spatial resolution of the SHEN- 
N (much larger than the landing site area) and uncertainties in the estimation of elemental 
composition. 

Figure 6. Smoothed (Gaussian filter with FWFIM = 120 km) map (in Robinson projection) 1° x 
1° of fast neutron flux (in counts per second) from the lunar surface measured by LEND SHEN- 
N detector. 



Figure 7. Average atomic mass of the sampled lunar soils versus fast neutron counting rate 
measured by LEND for the landing site regions. The best fit is shown by dashed line. 



[29] To compare LEND data and LPNS data, as was done for the thermal neutrons, we also have 
used latitude orbital phase profiles and estimated a correlation coefficient. The pixel-to-pixel 
correlation between LEND and LPNS maps is shown in Figure 8 . The linear correlation 
coefficient derived from this analysis is equal to ~0.6. It is significantly high value taking into 
account low counting statistics of SFIEN-N and the small regional variations of neutron flux 
everywhere except Fe enriched areas. The comparison of orbital phase profiles (OPPs) was also 
done and presented in Figure 9 . To create LEND and LPNS OPPs (we used LPNS map of fast 
neutrons provided in the special PDS data product based on the data reduction procedure 
described in Maurice et al. [20041) we used the technique described in section 3 . The LEND and 
LPNS orbital phase profiles follow each other. The small differences between them are covered 
by 1-er uncertainty of the LEND/LPNS normalized counting rates. Both OPPs show significant 
increases in counting rates (10% above average counting rate) at equatorial and moderate 
latitudes, corresponding to predominantly mare terrains, and lower counting rates in the 
highlands, confirming the correlation of fast neutron flux with average atomic mass distribution. 
At polar latitudes we do not see any substantial neutron suppression, as it is observed in 
epithermal neutrons (see next section). The comparison of epithermal and fast neutron fluxes 
distributions may be used for the verification of depth dependent distribution of Hydrogen at 
lunar poles. It may be done through detailed numerical modeling and use of a two layer model 
with an upper dry layer (with variable thickness) and a bottom layer with variable hydrogen 
concentration (probably in the fonn of dirty water ice). In analogy with Mars neutron data 


analysis, the HEND instrument integrated onboard the Mars Odyssey spacecraft [see, e.g., Litvak 
et al., 20061 shows that the thickness of the upper dry layer may be as large as 20 cm. 


Figure 8. The correlation between LEND SHEN-N fast neutron counting rate (y axis, counts per 
second) and LPNS counts from fast neutron detector (x axis, counts per second) using a 2° x 2° 
regular mapping grid. The linear best fit is shown by dashed line. 



Figure 9. Orbital phase curves (see section 4) showing regional latitude variations of fast neutron 
flux across the lunar surface. The x axis is the phase angle measured in degrees of arc from 0° up 
360°, where 90° corresponds to the south pole and 270° to the north pole. The y axis is the 
normalized counting rate of fast neutrons measured by different instruments. LEND data are 
shown with a solid curve and LPNS data with a dashed curve. 



5. Distribution of Epithermal Neutron Flux 

[30] In Figure 10 we show the map of counting rate detected by collimated sensors CSETNs (see 
Table 1) . Counting rates presented on this map are corrected for the detector efficiency changes 
(smooth saturation of counting rate after multiple switch offs/ons required by a correction of 
LRO orbit each 2 weeks) and for a long-tenn variations of GCR (variations of GCR flux as a 
function of solar cycle induce corresponding variations of Moon neutron flux). We also have 
taken into account a spacecraft background induced by the bombardment of the spacecraft body 
with particles of GCRs. The total value of background was estimated at the cruise stage far away 
from the Moon using the LEND measurements in the same configuration (detectors levels of 
high voltage, thresholds) as was set for the mapping phase. Then this value was adjusted for the 
LRO orbit (Moon partially shields spacecraft from the irradiation by GCRs leading to decreasing 
of the background) and subtracted it from the measurements. So, the result of the data processing 
presented in Figure 10 shows distribution of Moon neutrons only and logically can be divided 
into three different components. First one is the neutrons came into collimator FOV, second is 
the omnidirectional neutrons passed through collimator walls and finally third component 
presents omnidirectional neutrons backscattered in the spacecraft body. 

Figure 10. Smoothed (Gaussian filter with FWHM = 120 km) map (in Robinson projection) 1° x 
1° of counting rate (background subtracted) measured in LEND collimated detectors. 


[31] The map shows several distinctive regions with considerable deviations from the average 
counting rate (value averaged along all map pixels). First, polar regions display a depression of 
counting rate due to enhancement of hydrogen abundance at high latitudes (first revealed by 
LPNS). Second, at the nearside equatorial mare regions LEND CSETNs count significantly more 
neutrons than at the mid latitudes on the farside of the Moon. The corresponding quantitative 
latitude changes of CSETNs counting rate are presented in Figure 1 1 as a nonnalized orbital 
phase profile (same procedure as was used for thermal and fast neutrons in previous sections). At 
the poles it shows a 3-4% depression below average level of counting rate. On the contrary, at 
mare basins on the nearside of the Moon one can see 2-3% excess above average level of 
counting rate (phase angles 300°-360° in Figure 11) . We may conclude that spatial distribution 
of counting rate measured by LEND collimated detectors differs from the global maps created 
using data of omnidirectional neutron detectors. LPNS [see, e.g., Feldman et a A, 2001 ; Maurice 
et al., 2004] and LEND (SETN detector, see for description Table 1 and its normalized OPP in 
Figure 11) omnidirectional epithermal neutron detectors show less neutron counts at both poles 
(~4-5% of depression) and no significant excess of counting rate at the nearside mares (phase 
angles 300°-360° in Figure 11) . 

Figure 1 1 . Uncorrected orbital phase curves (see section 5) showing regional latitude variations 
of epithennal neutron flux across the lunar surface measured by LEND omnidirectional SETN 
(dash-dot line) and collimated CSETNs (solid line) detectors. Thex axis is the phase angle 
measured in degrees of arc from 0° up 360°, where 90° corresponds to the south pole and 270° to 
the north pole. The y axis is the normalized counting rate of epithermal neutrons. 


[32] In summary, we may conclude that the distribution of counting rates measured by LEND 
collimated detectors does not fully correlate with the distribution of epithermal neutron flux 
measured by omnidirectional neutron detectors. These distributions are similar at the poles 
showing evident suppression of counting rates but at midlatitudes on the nearside of the Moon, 
LEND collimated detectors count more neutrons than it is predicted from the omnidirectional 
detectors of epithermal neutrons (see Figure 11) . 

[33] Lawrence et al. [2011] have proposed that high energy epithermal and fast neutrons may 
pass through the collimator losing energy and shifting to the energy range where it may be 
detected by He tubes. Using relative variations of fast neutrons measured by LPNS they 
suggested that detection of high energy neutrons may add up to 2.25 counts per second into the 
total counting rate of the LEND collimated detectors. Using the same technique, Mitrofanov et 
al. [2011] have argued against this estimation and have shown that major “pollution” is 
contributed by fast neutrons and may be found as 1 . 1 counts per second. Mitrofanov et al. [2011] 
also have implemented more comprehensive approach comparing relative latitude variations 



(orbital phase profiles) of neutrons with different energies. It was done to understand how 
thermal, epithermal and fast neutrons latitudes variations should be mixed to determine latitude 
variations for counting rate observed in the collimated detectors. Following this approach it was 
again confirmed that penetration of fast neutrons describes the distribution of counting rates in 
the collimated detectors and amplitude of this component in the total counting rate is about 1 . 1 
counts per second. 

[34] In this paper we continue this analysis to derive a true global map of lunar epithermal 
neutron flux. 

[35] The basic idea of this method postulates that collimated and omnidirectional detectors 
(based on identical He tubes) may count different absolute numbers of neutrons (due to different 
FOVs) but relative variations of their counting rates should follow the same latitude profile 
shape. For example, if signal in a collimated detector changes by ~5% between latitude 1 and 
latitude2 (and the difference between these latitude belts is much larger then collimated FOV) we 
can predict the same relative change in the signal of the omnidirectional detector. It is obvious 
that the last statement is true only if both detectors count neutrons in the same energy range. The 
latitude variations for collimated and omnidirectional detectors can be presented as orbital phase 
profiles normalized to the average counting rate measured by each detector. Applying OPP for 
the epithermal LEND data sets (CSETNs and SETN) we could expect that OPPcsetns(^) was 
equal to the OPPsetn(^)- In reality, as we already established (see Figure 1 1) , their shapes are 
quite different. 

[36] To explain this difference we propose that the total counting rates in CSETNs and SETN 
may also include contributions from thermal, high-energy epithermal and fast neutrons. The 
CSETNs are well protected from thermal neutrons (by collimator body, additional cadmium 
and 10 B sheets on the detectors tops) but it may be sensitive for the high-energy epithermal and 
fast neutrons moderated in the collimator walls or backscattered in the spacecraft body. 
Comparing latitude variations of CSETNs and SETN (presented in Figure 11) one may state that 
“non-epithermal” component in the CSETNs total counting rate both increase neutron flux at the 
poles (because poles depressions of CSETNs curve are significantly smaller than SETN ones) 
and at the equatorial latitudes at the nearside mares (CSETNs curve lies above SETN one at the 
range of phase angles 300°-360°). The best model of such “non-epithermal” component is the 
latitude variations of the fast neutrons measured by LEND SHEN-N detector (see section 4) . 
According to Figure 9 the OPP for SHEN-N detector has the significant increase of counting rate 
at the midlatitudes (phase angles 300°-360°) and negligible depression at the poles. 

[37] On the other side, SETN (like LPNS epithermal detector) is not sensitive to the high energy 
and fast neutrons [see, e.g., Lawrence et al., 201 1] but may detect some small portion of the 
thermal neutrons leaking through the cadmium enclosure of SETN. 

[38] Concluding this discussion we may say the following: To attain the true counting rate of 
epithennal neutrons in the collimated detectors we need to subtract the unknown fraction of 
counting rate from the detection of fast neutrons that propagate through the collimator walls or 
backscattered in the spacecraft body. The latitude behavior of this “background” component is 
described by OPPshen-n- To get the true counting rate of epithermal neutrons in the 


omnidirectional detector SETN we possibly need to subtract the u nkn own fraction of the 
counting rate deposited by thermal neutrons. The latitude behavior of this component can be 
described by OPPooppier described in the section 3 (see Figure 4) . 

[39] Mathematically it may be presented as two equations: 

. = ( 1 - a) x OPP?£r' + a x OPP^. N (3) 

. OPPgg, = (1 - 0) x OPPST 1 +0* OPPgZZ (4) 

where OPP lolal esETNs and OPP Total sETN are relative latitude variations measured by CSETNs and 
SETN (presented in Figure 1 1) ; a and fi are free parameters; OPPshen-n- normalized counting 
rate derived from LEND fast neutron detector SHEN-N (see section 4 and Figure 9) and 
OP P Dopplers is a difference of counting rate measured by the Doppler detectors presenting global 
variations of thermal neutrons (see section 3 and Figure 4) . OPP Epitherma ’csETNs and 
OPP Ep ithermai SETN are normalized true epithermal counting rates in CSETNs and SETN. 
Postulating that they must be equal each to the other, we may solve the system of equations and 
find the best fit values of a and [:> for the whole range of phase angles, >p. Quality of this fit is 
verified by a Pearson criteria test based on a least squares analysis. We have used expression 

(5) 


where A is a OPP Epithermal csETNs expressed from equation (3) as a function of a, B is 
OPP E pithermai sETN cx p rcsscc [ f rom equation (4) as a function of /?, S- is a total uncertainty derived 
from the statistic of counts in CSETNs, SETN, Dopplers and SHEN-N detectors, {^i} - set of 
phase angles. Varying a and [i we may minimize S 2 . According Pearson criteria if difference A-B 
agrees within statistical uncertainties then S“ shall follow distribution/ with N degree of 
freedom where N is equal to number of summands in sum (5) reduced by number of free 
parameters \ Press et al . , 19921 . 

[40] We have discovered that it is possible to find such values of parameters a and /?, which will 
be accepted by Pearson criteria. It confirms that suggested division on the components was 
reasonable. The results of this analysis are given in Figure 12 where we have shown corrected 
OPP for the collimated detectors (corrected for high energy/fast neutrons) and OPP for the 
omnidirectional detector SETN(corrected for the thermal neutrons). This analysis leads us to an 
estimation of the best fit values of fractions of “true epithermal neutrons” and “non-epithermal 
neutrons” measured by the collimated and omnidirectional detectors. So, the resulting counting 
rate in the collimated detectors produced only by epithermal neutrons in the FOV is equal 
to ~ 1 .7 counts per second, in good agreement with our previous estimates [ Mitrofanov et al . , 
2010a , 2010b , 2011 ; R. S. Sagdeev et al., Use of Apollo 17 epoch neutron spectrum as a 
benchmark for the Collimated 1 Sensors of the Lunar Exploration Neutron Detector, submitted to 
Journal of Geophysical Research, 2012]. The counting rate of epithermal neutrons in SETN is 
about 9.8 counts per seconds. Farther we will use this signal to present a global map of 
epithermal neutrons. 


Figure 12. Corrected orbital phase curves (see section 5) showing regional latitude variations of 
epithennal neutron flux across the lunar surface measured by LEND omnidirectional SETN 
(dash-dot line) and collimated CSETNs (solid line) detectors. Thex axis is the phase angle 
measured in degrees of arc from 0° up 360°, where 90° corresponds to the south pole and 270° to 
the north pole. The y axis is the normalized counting rate of epithermal neutrons. 



[41] Figure 13 shows OPP for the SETN signal corrected for the “non-epithermal” component 
(same as in Figure 12) and OPP created from the LPNS PDS maps data sets (special data 
product), but with the following difference: the epithennal neutron counting rate is reduced by 
£*thermal neutrons, where e is equal to 0.057 as recommended in Feldman et al. [200 1] , 
Comparing these OPPs one may conclude that the shape and amplitude of neutron extended 
suppression measured by LEND around the lunar poles agrees with the LPNS data. But one can 
see that at equatorial highlands on the farside of the Moon, the LEND average counting rate is 
higher than that given by LPNS. At moderate latitudes in the north, LPNS predicts a higher 
intensity of epithennal neutron flux compared to LEND data. It can be seen in more detail in 
Figure 14 where we have shown a map of relative differences between LEND and LPNS 
epithermal neutron fluxes expressed in the statistical deviations (same approach as used in the 
section 3 , equation (1)) . The yellow and red areas on this map correspond to maximal 
discrepancies between LEND and LPNS data. The largest discrepancies between LEND and 
LPNS are seen in the northern hemisphere at Mare Frigoris where LPNS has significantly higher 
values of epithermal flux than LEND and at equatorial latitudes [ — 30S, 30N] in the highlands on 
farside (with longitudes from -180° up to -120° and from 120° up to 180°) of the Moon where 
in contrast LEND is counting more epithennal neutrons than LPNS. In Figure 15 we also have 
presented LEND/LPNS pixel-to-pixel correlation. It shows fairly well linear correlation 
coefficient is about 0.73. It is easily seen (in Figure 14) that the best agreement between LEND 
and LPNS data is in the southern hemisphere poleward of ~50 S for all longitudes. In the north 
there is a good agreement near polar latitudes (the same effect of extended polar suppression). 
The significant correlation is also found poleward of 60 N within [-180°, -120°] and [60°, 180°] 
longitude sectors. The same pattern of distribution of low values of epithermal flux is seen at 
moderate and near equatorial latitudes in the vicinity of Mare Imbrium basin. 

Figure 13. Corrected orbital phase curves (see section 5) showing regional latitude variations of 
epithermal neutron flux across the lunar surface measured by LEND omnidirectional SETN 
(dash-dot line) and LPNS (solid line) detectors. Thex axis is the phase angle measured in degrees 
of arc from 0° up 360°, where 90° corresponds to the south pole and 270° to the north pole. The 
y axis is the normalized counting rate of epithermal neutrons. 
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Figure 14. The map of comparison (number of statistical deviations) between LEND epithermal 
neutron counting rate (SETN detector) and epithermal neutron counts from LPNS detector. For 
details, see section 5 and equation (1) . 



Figure 15. The correlation between LEND corrected epithermal neutron counting rate (SETN, y 
axis, counts per second) and LPNS counts from epithermal neutron detector (xaxis, counts per 
second) using a 2° x 2° regular mapping grid. The counting rate in SETN was corrected for the 
detector efficiency changes, long-term variations of GCRs. The counting rate was also reduced 
by the GCR spacecraft background (defined at the cruise) and by contribution of thermal neutron 
counting rate derived from equations (3) and (4). The linear best fit is shown by dashed line. 



[42] The resulting corrected global and polar maps of LEND epithermal counting rate (as 
measured by omnidirectional sensor with spatial resolution, —75 km and reasonably good 
counting statistics, average relative error <2%) are shown in Figures 16 and 17. They are done 
with a 1° x 1° grid and smoothed with a Gaussian filter with FWFIM equal to 60 km (compared 
with omnidirectional spatial resolution). Several large crater basins both in the farside and 
nearside show increasing intensity of epithermal neutron flux in comparison with neighboring 
areas. For example, one can see a substantial increase of epithermal neutron flux in the southern 
hemisphere in such crater basins as Clavius, Maginus, Longomontanus and Tycho (marked in 
Figure 16) . At the equatorial latitudes within [-30S, 30N] high epithermal flux is observed at 
Mare Orientale, Mare Fecunditatis, Mare Nectaris, Mare Crisium and Mare Moscovience. 

Figure 16. Smoothed (Gaussian filter with FWFIM = 60 km) map (in Robinson projection) 1° x 
1° of epithermal neutron counting rate (background subtracted, thermal neutron component 
removed, counts per second) measured in LEND SETN detector. The counting rate in SETN was 
corrected for the detector efficiency changes, long-term variations of GCRs. The counting rate 
was also reduced by the GCR spacecraft background (defined at the cruise) and by contribution 
of thermal neutron counting rate derived from equations (3) and (4) . 



Figure 17. (top) Smoothed (Gaussian filter with FWFIM = 60 km) map (in stereographic 
projection) 1° x 1° poleward of 70 S of epithermal neutron counting rate (background subtracted, 
thermal neutron component removed, counts per second) measured in LEND SETN detector. 

The counting rate in SETN was corrected for the detector efficiency changes, long-term 
variations of GCRs. The counting rate was also reduced by the GCR spacecraft background 


(defined at the cruise) and by contribution of thennal neutron counting rate derived from 
equations (3) and (4). (bottom) Smoothed (Gaussian filter with FWHM = 60 km) map (in 
stereographic projection) 1° x 1° poleward of 70 N of epithennal neutron counting rate 
(background subtracted, thermal neutron component removed, counts per second) measured in 
LEND SETN detector. The counting rate in SETN was corrected for the detector efficiency 
changes, long-term variations of GCRs. The counting rate was also reduced by the GCR 
spacecraft background (defined at the cruise) and by contribution of thermal neutron counting 
rate derived from equations (3) and (4). 



[43] Polar maps from omnidirectional detectors cannot help us to distinguish correlations of 
neutron suppression with distribution of polar permanently shadow craters due to poor spatial 
resolution. But polar maps of epithermal neutron flux measured in LEND collimated sensors can 
be created on a much finer grid. These maps may be used to search for local areas with counting 
rates showing significant neutron suppression or excess. The detailed discussion concerning such 
an analysis is given in I. G. Mitrofanov et al. (Testing polar spots of water-rich permafrost on the 
Moon: LEND observations onboard LRO, submitted to Journal of Geophysical Research, 2012); 
W. V. Boynton et al. (High spatial resolution of epithermal neutron emission from the lunar 
poles: Constraints on hydrogen mobility, submitted to Journal of Geophysical Research, 2012; 
and Sanin et al. [20121 (distribution of epithennal neutron flux inside and outside Polar 
Shadowed Regions). 

6. Summary 

[44] We have presented the latest results of global mapping of lunar neutron flux measured by 
the LEND instrument onboard the LRO mission. Three neutron energy ranges, thermal 
(<0.015 eV), epithermal (0.4-500 eV) and fast (>0.5 MeV) have been studied in our analysis. 
Thermal neutrons are sensitive to specific elements in the lunar subsurface with high absorption 
cross sections (like Fe). Hydrogen is a dominant factor in determining the variations of 
epithermal neutron flux. The changes of fast neutron flux mainly follow the average atomic mass 
of lunar regolith. Global maps of derived (corrected for systematic effects) counting rates have 
been separately created for thermal, epithermal and fast neutron energy ranges. The regional 
variations observed on these maps are reasonably well correlated on a global scale with previous 
(LPNS) results, but show differences on a local scale. 

[45] The largest regional variations (up to a factor of 3-4) are seen in the LEND map of thermal 
neutrons counts between the Fe-rich mare basaltic terrains on the nearside and Fe-poor highlands 
on the farside of the Moon (see Figure 2) . It is similar to the LPNS results [see, e.g., Lawrence et 
al, 20021 . 


[46] The variation of the epithennal neutron flux across the lunar surface (SETN map, Figure 16) 
is much smaller (~ 10%) and visible in both LEND and LPNS data as the extended neutron 
suppression poleward of 70S and 70N. At moderate and equatorial latitudes (where observed 
regional variations are much smaller) the differences between LEND and LPNS maps are most 
noticeable. It may be a result of counting statistics or different approaches in data reduction 
procedures. 

[47] The fast neutron flux changes by 25% across the lunar surface from maria to the highlands, 
with the maximum of fast neutron flux observed in mare terrains. This result is consistent with 
the LPNS argument that mare basalts are rich in Fe (producing more fast neutrons in comparison 
with Al-rich soils). The latitude band profile of fast neutron flux does not show a significant 
polar extended neutron suppression effect as observed in the epithermal neutron range. One may 
distinguish only small areas around poles with neutron suppression less than 1%. By comparison 
to detection of water ice at Mars polar latitudes by HEND [see, Q.&. Xitvak et al., 2006] , this may 
be interpreted as the hydrogen distribution at the lunar poles, on average, is depth dependent with 
higher weight fraction of H beneath a relatively hydrogen poor regolith. 
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